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We present limits on the ejection of old-population HVS from a sample of 



over 290,000 stars selected from the SDSS. We derive the speed at the solar circle 
from the measured positions and radial velocities by assuming a radial orbit and 
adopting a simple isothermal model of the Galactic halo, which enables us to 
identify candidate bound and unbound ejectees. We find 4 candidate bound 
F-stars from this sample, all with negative Galactocentric radial velocity (i.e., 
returning toward the GC). We additionally find 2 candidate unbound stars (one F 
and one G), however, existing proper motion measurements make these unlikely 
to be emerging from the GC. These data place an upper limit on the rate of 
ejection of old-population stars from the GC of ~ 45Myr _1 . Comparing to the 
rate for more massive B-star ejectees of ~ 0.5 Myr -1 , our limit on the rate of 
ejection of old-population HVS shows that the mass function at the GC is not 

1 1 

bottom-heavy and is consistent with being normal. Future targeted surveys of 
^ ! old-population HVS could determine if it is indeed top-heavy. 
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1. Introduction 

The search for hypervelocity stars (HVS) in the stellar halo was original ly propose d 

S3' 



in order to indirectly probe the Galactic Center (GC) at optical wavelengths faills 1988 ). 



These stars, having been ejected from the GC via binary-exchange collisions with a putative 
supermassive black hole at speeds in excess of Galactic escape would provide the dynamical 
"smoking gun" evidence for such an object at a time when direct detection seemed remote. 
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Despite the tremendous advances in infrared instrumentation in the intervening two decades 
(and the associated confirmation of the supermassive black hole at the GC), it remains 
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Galaxy (e.g. 

of these young stars is a major mystery since one would naively expect that the conditions 
near a supermassive black hole are sufficiently hostile to prevent in situ star formation of any 
kind. A key question is whether or not these stars represent the "tip of the iceberg" in an 
otherwise normal star forming region or whether star formation is fundamentally different at 
the GC. Because HVS are relatively rare, it was previously not feasible to amass sufficiently 
large numbers of such stars to go beyond the "smoking-gun" prediction that there should 
be young, blue stars in the halo. In the era of deep, wide-field spectroscopic surveys, it is 
now conceivable to use Hills' original idea to exploit HVS in the Galactic halo to obtain, 
among other things, detailed information about star formation at the GC, which is otherwise 
obscured from direct vie w, to at least a factor of 2 .5 lower in mass than the young population 
currently being probed (IBrown et al.ll2005l 120081 ). 



The Sloan Digitial Sky Survey (SDSS) provides an exceptional database to search for 
HVS that would otherwise be missed from targ eted surveys. The SDSS has obtain e d ove r 
300,000 stellar spectra since it began operations (jYork et al.ll2000l ; iGunn et al. Ill998l . 120061 ). 
The data are reduced by an automated pipeline that produces reliable spectral classifications 
and radial velocities (RVs). The SEGUE Stellar Parameter Pipeline (hereafter, SSPP) fur- 
ther processes the wavelength- an d flux-calibrated spect ra generated by the standard SDSS 
spectroscopic reduction pipeline (jStoughton et al.l 120021 ). obtains equivalent widths and/or 
line indices for 82 atomic or molecular absorption li nes, and estimate s T e ff, log q, and [Fe/Hl 
throu gh the application of a number of approaches kee et. riWJfl d -Prieto et al.l 



20081 ). With this enormous database of stellar spectra it is already possible to detect HVS 
among a broad class of stars. In particular, the size of the SDSS stellar library makes SDSS 
sensitive to F and G type H VS which require a subst antially different search strategy than 
their O and B counterparts (iKollmeier &; Gouldl 120071 ) . In contrast to the photometric sur- 
vey, the SDSS RV catalog is not complete in any dimension. In this work, we examine a 
large sample of high- velocity stars in order to place limits on the ejection of F/G stars from 
the GC. 

In § [2] we will summarize the stellar spectra and the sample selection for stars that will 
be used for this analysis. In § [3] we present limits on the ejecton of old-population HVS. 
We compare this limit with young-population HVS in § |H We present our discussion and 
conclusions in § 
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2. The Sample 

The sample used here is part of an ongoing study to find high-velocity stars within the 
SDSS (Knapp et al. 2009, in preparation). All RVs for stars aquired by the SDSS as of 15 
January 2007 were extracted from the SDSS database. Selected objects were required to have 
RV errors not exceeding 100 kms" 1 and required to be free of pipeline flags indicating likely 
problems in the automated velocity determination. This yielded a sample of 291,111 objects. 
Velocities were converted to Galactocentric velocity assuming that the local rotation speed 
is 220 kms -1 and that the Sun moves at (+9, +12, +7) kms" 1 relative to the local standard 
of rest in the direction of the GC, Galactic rotation and the north Galactic pole respectively. 
All objects with Galactocentric RVs in excess of \Vq\ > 350 kms" 1 were examined to ensure 
the sample was free of catastrophic velocity errors. This process yielded a total sample of 33 
objects with Galactocentric RVs exceeding ±400 kms" 1 , which we analyze in detail below. 



3. F and G Star Ejectees 



Kollmeier fc Gouldl (120071 ) argued that large RV samples, such as SDSS, could probe an 
as yet undiscovered old population of stars ejected from the GC and would be most sensitive 
to stars near the main-sequence turnoff within this population. Moreover, they argued that 
such surveys would be more sensitive to the bound rather than unbound members of this 
population, simply because the bound stars accumulate over the lifetime of the Galaxy and 
can be seen still orbiting, while the unbound stars can only be viewed during their exit. In 
this section, we analyze the sensitivity of SDSS to both classes of ejectees, and we tabulate 
the candidates derived from this database. 



3.1. Selection Criteria for Bound Ejectees 

We begin by asking what the sensitivity of the SDSS survey is to bound turnoff stars 
and whether any such candidates are in the sample. Of course, this requires that we establish 
selection criteria that remove the vast majority of contaminants but still retain significant 
sensitivity to bound ejectees. We adopt the following two criteria: 

V r ,G > V cnt = 400 kms" 1 , V -circle < V esc , (1) 

where t> r ,G is the observed RV converted to the Galactic frame, v ©-circle is the velocity that 
the star has when it crosses the solar circle, and t> esc is the Galactic escape velocity, again 
measured at the solar circle. 
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The first criterion is purely observational (apart from the implied assumption of solar 
motion) and is self-contained. The threshold i> cut is established empirically from the observed 
velocity distribution (Knapp et al. 2009, in prep). The second criterion requires that we 
specify both a model of the Galaxy (to determine the value of v esc ) and a model of the 
photometric properties of the target population (ejectees), so that we can estimate w _ C i rc i e 
from its observed tv,G and from the observed fluxes in several bands. 

For the Galactic model, we adopt an isothermal sphere characterized by a rotation 
speed v c = 220 kms -1 , truncated at a Galactocentric radius Rf,. For this model = 
Ro exp[(t> esc /t> c ) 2 /2 — 1]. We will adopt v esc = 550 kms -1 , which implies Rb = 8.37-Ro, as a 
default, but will also consider other values. We will consider two classes of stars,' "F stars", 
defined observationally as 0.3 < (g— i) < 0.5, and "G stars", defined as 0.5 < (g— i) < 0.75. 
We adopt for these absolute magnitudes M g = 3.5 and M g = 5.5, respectively. It is important 
to keep in mind that these assumptions apply only to the putative ejectee population: the 
photometric properties of non-ejectee stars, which obviously dominate the SDSS sample 
overall, are completely irrelevant. 

From the observed dereddened magnitude go of the star, we can then infer its distance, 
r = \^9o-M a )/r>-2 kp C ^ an( j j ience (f rom its position on the sky), infer the angle between 
the Sun and the GC, as seen from the star. Then, since the ejectee is assumed to be on a 
radial orbit, its full 3-D velocity is v — |tv,G sec 0|, and hence its velocity at the solar circle 
(assuming, as is always the case for F stars in SDSS, that R < R b ) is 



It is immediately obvious that this equation places strong constraints on the observational 
parameter space in which the search can be conducted. For instance, if the star is fairly 
bright, so that r «C Ro, then cos0 > v cut /v esc . Since, for this geometry, cos0 ~ — cos /cos b, 
this constraint eliminates the roughly f cu t/^esc = 73% of the sky that is not sufficiently close 
to the Galactocentric (or anticentric) directions. At the opposite extreme, a bound star at 
R > Ro exp[(i7 csc /t> cut ) 2 /2] ~ 2.57 Ro cannot be probed in any direction because its velocity 
(and hence RV) will fall below the threshold. At intermediate distances, these two effects 
combine with different relative weights. 



Of the stars in the Knapp et al. (2009) sample, four F stars (and no G stars) survive 
as candidate bound ejectees (see Table [1]). It is a curious fact that all four are coming 
toward us, however there is a 1/8 random probability that all four would be going in the 




(2) 



3.2. Candidate Bound Ejectees 
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same direction, which is too high to warrant further analysis. It is also interesting that these 
stars are all measured to be metal poor. These candidates could either be members of a 
metal-poor population ejected from the GC or simply halo stars. Obtaining proper motion 
measurements for the candidates would discriminate between the two scenarios. 



3.3. General Formulae for Bound-Ejectee Sensitivity 

Initially, consider a spectroscopic survey that obtains RVs of all stars in a small angular 
area Q over a narrow magnitude range go±Ago/2, and over a specified narrow range of colors. 
Consider a star in this sample that is ejected from the GC, with current Galactocentric RV 
v r , and assume that it has an absolute magnitude (estimated from its color) M g . The star 
then has distance from us r = 10^ ()_M9 ^ 5-2 kpc, so the volume of such a survey is 

AV = ^% fir 3 . (3) 

5 

As mentioned above it has a 3-D Galactocentric velocity v — v r sec 0, where is the angle 
between us and the GC, as seen from the star. 

As seen from the GC, the probed volume covers an area AA and has a thickness Aw. Of 
course, AV = AAAw. A fraction AA/AirR 2 of all ejected stars will pass through this volume 
at some time (assuming isotropic ejection), where R is the star's Galactocentric distance, 
and it will spend a fraction 2 Aw / vP of its time in this volume, where P is the orbital period. 
Hence, observations of this volume will probe a fraction 

AA 2Aw _ In 10 Ag ttr 3 
^ ~ AnR 2 vP ~ 2.5 AttR 2 vP ^ ' 

of all ejected stars. Now consider that we do not observe all stars in this volume, but only 
one. Clearly, the fraction of ejected stars that we probe falls by a factor iV vo i, where iV vo i is 
the number of stars that satify our detection criteria. That is, N vo \ = nflAg , where n is 
the number density of stars per sterradian, per magnitude, and satisfying the color criteria. 
Hence, the fraction of all ejected stars probed by this single observation of the iih star is 

In 10 

Jl 2.5 A-KRlvPiTii y ' 

To determine the period, we apply our adopted Galactic model, first finding that the 
apocenter of the orbit R m will be 

R 

R m = i?*, (R*<R b ); R m = 1 + HRb/R) b _ {v/Vc)2/2 (R*>R b ) (6) 



- 6- 



where 



i?* = i?exp 



{v/v c f 



The velocity v' at any point in the orbit R' will be given by 

2\v r J R r v - ; ' 2\vJ R' Rrn 



(7) 



(8) 



Integrating the equation dt = dR'/v'(R') over one full period yields 

-R 



P = v^^erfc Mln(R*/R b )] + v^— [v 7 !^ + 
where x = R b / R m . 



arccos \/x 



x 



(9) 

(R* > R b ) (10) 



The total sensitivity of the survey is then simply the sum / tot = £\ /j over all stars. If 
there are iV tot bound stars orbiting, then the expected number of detections will simply be 



-ex P — /tot-^tot- It will prove useful later on to express the same result as 

^thresh = T MW ^ ft, 



(11) 



where tmw — 10 Gyr is the lifetime of the Milky Way. In this formulation, A^ exp = r/r thresh , 
where T is the mean bound-ejection rate averaged over the lifetime of the Galaxy. 

To understand the basic properties of f iy let us initially restrict attention to the first 
case, i?* < Rb- Then 



-3/ 2 /^ 3 exp(-^72) 
5 v "" / U 



/ 4 = ^Vr^(-f 



zn.: 



V 



(12) 



If we consider some typical values for relatively faint turnoff stars in SDSS (g ~ 18.5), 
r ~ lOkpc, R ~ 16kpc, z — 2, and n = 100 deg -2 , then ~ 1.5 x 10~ 9 . Hence, by 
obtaining RVs for about 10 4 stars, one could probe a bound population of F stars ejected 
from the GC provided it had at least 10 5 members. This corresponds to r throsh = lOMyr -1 . 



3.4. Bound-Ejectee Sensitivity of SDSS 



To determine the sensitivity of the SDSS spectroscopic survey, we first evaluate n» for 
our two classes of stars as a function of magnitude and position on the sky using the SDSS 
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photometric catalog. This was done using 20 individual patches of sky, each with area 1 
square degree, and interpolating these values to obtain n.j for other parts of the sky. We 
then sum equation (jSJ) over all stars in the SDSS spectroscopic catalog that meet our color 
criteria. Figure [T] shows the result for several cases. The three bold curves (blue) are for F 
stars under the assumption of t> esc = 550 kms -1 and for three different thresholds, t> cut = 400, 
410, and 420. The solid (red) curve is for F stars with t> cut = 400 kms -1 , t> esc = 520 kms -1 , 
while the dashed (green) curve is for G stars with t> cut = 400 kms -1 , t> esc = 550 kms -1 . The 
two other curves will be explained later. 

Several features are apparent from this figure. First, sensitivity is peaked quite close to 
the escape velocity. This may seem surprising in view of equation (fT2"]) . which exponentially 
declines with the square of the velocity for R* < Rb, and then even more steeply. However, 
at low v, distant (i.e., faint) stars do not contribute because their declining speeds at large 
R bring them below the first selection cut in equation ([1]). Most distant stars gain by r 3 
in equation (fT2"j) . at least for r < R . Second, sensitivity falls quite steeply with increasing 
^cutj approximately a factor 2 for each 10 kms -1 . This is basically a consequence of the 
same volume effect just analyzed. The same effect also accounts for the dramatic falloff in 
sensitivity if the escape velocity proves lower than our default f esc = 550 kms -1 . Finally, 
senstivity to bound G-star ejectees is about 1/3 that of F stars. 

Let us initially adopt the extremely simple model that the number of bound F-star 
ejectees is uniform as a function of f Q _ circ ie. Under this assumption, we can make one test of 
the plausibility of the candidate population by ploting their cumulative distribution, which 
is shown in Figure [2j The match is extremely good. 

Let us then tentatively assume that the candidates are really bound ejectees. Inte- 
grating under the entire (upper-blue) curve in Figure [T] yields J dvT^.^ = 4.0 Myr- km s -1 . 
However, from Figure [T] and especially Figure [2J it is clear that about 75% of this sensitiv- 
ity comes from the roughtly 35 kms -1 interval in which the four candidates are detected. 
Therefore, if these candidates are accepted as real, we can infer an ejection rate 

4 35 kms -1 

45 Myr -1 (505 kms -1 < v _ circ i c < 540 kms -1 ) 



<r -1 resh ) 0.75x4Myr -1 -kms -1 

(13) 

Of course, if some or all of the candidates are later vetted and found wanting, the smaller 
number of detections would then yield an analagous upper limit. 



Since the G-star sensitivity is smaller and we find no candidates, the G stars are con- 
sistent with a similar rate, but with weaker statistics. 
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3.5. Unbound Ejectees 

We select candidate unbound ejected stars with exactly the criterion as bound stars 
(eq. pQ), except reversing the sign on the second condition, %- c ircie > v esc , and of course 
demanding that they be exiting rather than entering the Galaxy. Then, following a very 
similar derivation to the one given for equations (jSJ) and (|lip . we immediately derive their 
analogs for the unbound case, 

r- 1 =r f . t - ln 10 rf (u) 

1 thresh 2> 5 A7T R2 vn .- ^ 

i 

This quantity is shown in Figured] as a function of t> Q _ C i rc i e for F stars (bold dashed, magenta) 
and G stars (dotted, cyan) for the case f esc = 550 kms -1 , t> cut = 400 kins' 1 . Nothing 
prevents us from evaluating equation (TT4"|) even below t> esc , and we display the extensions of 
these curves in Figure [TJ Their physical meaning is that it is possible to detect bound ejectees 
during their first orbit, and for orbital times of order or longer than tmw, this is actually a 
more accurate representation of the detection rate than averaging over a full orbits. Thus, 
the bound rate is really the maxium of the bound curve and extension of the unbound curve. 

There is one F star candidate that survives our selection criteria and one G star. The 
proper motions for both of these stars are inconsistent with a Galactocentric origin as we 
now show. Were it ejected from the GC, based on its current position, a star should exhibit 
a proper motion in the direction 

tan(0) = tan(/)sin(6) (15) 
where 9 is measured from G alactic north through east. The F st ar candidate has a measured 



proper motion from SDSS ( iPier et al. 1120031 : iMunn et al.ll2004j ) of 34masyr 1 at 55° (with 



small errors) — significantly in conflict with the predicted value 9 = 330°. 

Unlike the F star, the G star candidate cannot be discarded based on the its proper- 
motion direction alone. The measured proper motion of 6masyr _1 is sufficiently small 
that the measured direction of 183° is poorly determined given the errors. We therefore 
compare the predicted and observed vector proper motions, which (unlike the case of the F 
star) require a distance estimate. From its g magnitude we estimate a distance of 7.8 kpc, 
and it is therefore expected to be moving at 14masyr _1 at 244°. The magnitude of the 
vector difference of the expected and observed p roper motions is l^ masyr -1 , which is to be 



compared to the typical errors of ~ 4 mas yr 1 (IMunn et al.l 120041 1 , making it unlikely that 
this star is an unbound ejectee. 

Since Thresh is essentially constant in the unbound regime at (0.05 Myr)" 1 for F stars 
and (0.01 Myr) -1 for G stars, this lack of detections places an upper limit on the ejection of 
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unbound F and G stars of 

rebound < 60 Myr-\ T^ hound < 300 Myr^ 1 , (16) 

at 95% confidence. 

Equation (flBl) is broadly consistent with equation (fl3|) . 

4. Comparison with B stars 

The Brown et al. (2007) survey in essence covered an entire volume of Galaxy rather 
than a scattering of targets within the Galaxy, like SDSS. Thus, to calculate Thresh one 
should perform a volume integral. An argument similar to the one given in § 13.31 leads to 
the equation 

! f dsinbdl r max drr 2 

thrCSh = Lea-covered X min I W P " ^1 M^o)] ^ ^ 

where we have for convenience assumed the form of the Galactic potential within < Rb. We 
esitmate this expression for the Brown et al. (2007) survey by noting that it covered an area 
AQ = 5000 deg 2 and that most of the sensitivity was in regions where R ~ r, by adopting an 
average value for the velocity in the denominator of (v) = 600 kms -1 , and approximating 
Ar = r max - r min = 40kpc. We then obtain T~^ resh = (AQ / 4n) (Ar / (v)) = 8Myr. Hence the 
four detections imply r u ? nbound = 4r th r es h = 0.5 Myr -1 . 

Thus, the upper limit given by equation (TTBTl shows that F-star ejectees are no more 
than 100 times more common than B-star ejectees. Since F stars themselves are about 100 
times more common than B stars in the Galaxy as a whole, this result is already of some 
interest. Nevertheless, it would certainly be better to achieve at least a few times better 
sensitivity. 



5. Discussion 

We have presented the first limit on old-population ejectees from the Galactic Center 
based on data obtained with the SDSS. We have shown that the ejection of F/G stars is 
no more than 100 times more common than O/B star ejectees. We note that this limit is 
obtained simply from analyzing the stellar spectra obtained by SDSS as part of its main 
projects, not via a survey designed to preferentially select these objects. Already, this limit 
suggests that the mass function at the Galactic center is not weighted toward low mass star 
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formation — a fact impossible to determine by direct imaging of the GC alone. Because 
these candidates are all found to be metal poor (see Tabled]), they may be part of the metal- 
poor halo in which case, this limit is even stronger. Proper motion measurements would 
determine whether they are indeed ejectees from the GC and thereby make this limit more 
stringent. Ongoing surveys, such as the second epoch of the Sloan Extension for Galactic 
Understanding and Exploration (SEGUE-II), will determine conclusively if star- formation 
at the GC is indeed weighted heavily toward young, massive stars. 
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Table I. Candidate Ejectees Satisfying Velocity Criteria 



Name 


Vq (km/sec) 


Type 


9Q 


(g - i)o 


Toff (K) 


log(g) 


[Fe/H] 


Plate/Fiber/MJD 


SDSSJ061118.63+642618.5 


-412.3±5.5 


F 


18.39 


0.45 


6143 


4.08 


-1.27 


2299/488/53711 


SDSSJ074557.31+181246.7 


-409.0±5.9 


F 


18.96 


0.34 


6235 


3.51 


-1.73 


2074/113/53437 


SDSSJ224052.56+011332.1 


-407.1±16.3 


F 


19.30 


0.49 


6168 


4.02 


-1.37 


1101/561/52621 


SDSSJ211321. 02+103456.7 


-404.4±17.2 


F 


20.07 


0.35 


6069 


4.10 


-1.91 


1891/423/53238 


SDSSJ111107.85+585357.2 


+425.4±3.1 


F 


16.16 


0.36 


6514 


3.58 


-1.77 


950/554/52378 


SDSSJ224740.09-004451.6 


+401.8±13.1 


G 


19.97 


0.72 


5518 


4.17 


-1.22 


1901/7/53261 



450 



500 550 
v at Sun (km/s) 



600 



Fig. 1. — Sensitiviy of SDSS to old-star ejectees as a function of type and Galactic model. 
Bold, blue curves show the sensitivity for bound F stars assuming t> esc = 550 kms -1 and 
for three different thresholds, t> cu t = 400, 410, and 420 from bottom to top. Solid red curve 
shows the case for bound F stars with v cut = 400 kms -1 , t> esc = 520 kms -1 , while the dashed 
(green) curve is for bound G stars with t> cut = 400 kms -1 , v csc = 550 kms -1 . The sensitivity 
to unbound F stars (bold dashed, magenta) and unbound G stars (dotted, cyan) for the case 
f esc = 550 kms -1 , f cut = 400 kms -1 are also shown. 
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Fig. 2. — Cumulative distribution of velocities of candidate bound F-star HVS population 
(red histogram). A simple model in which the number of bound F-star ejectees is a uniform 
function of v ©-circle is shown in black for comparison. 



